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Introduction 

Extending the NASA Ames Mars General Circulation Model (MGCM) upper boundary will expand our 
understanding of the connection between the lower and upper atmosphere of Mars through the middle 
atmosphere. 

The extension's main requirements is incorporation of Non-local thermodynamic equilibrium (NLTE) heating 
(visible) and cooling (infrared). 

NLTE occurs when energy is exchanged more rapidly with the radiation field (or other energy sources) rather 
than collisions with other molecules. 

Without NLTE, above ~80km/~60km in Mars' atmosphere the IR/visible heating rates are overestimated. 
Currently NLTE has been applied successfully into the ID RT code and is in progress for the 3D application. 


ID Radiative Transfer Code 

RT code is based on two-stream approximation; quadrature in the visible, hemispheric mean in the IR 
Calculates fluxes at layer boundaries and at the surface; heating rates are calculated from these fluxes 
Includes C0 2 and H 2 0 gas opacity; aerosol optical properties for dust and water ice clouds 
Correlated-k gas opacities; calculated off-line on pressure and temperature grid 
Pressure ranges: 10 -6 to 10 +4 mbar (pressure at each grid point lOx larger than previous) 

Temperature ranges: 50 to 350 K, AT = 50 

12 spectral intervals (7 visible, 8 IR) with 17 Gauss points per interval 
IR spectral range: 1000 4.5 pm; Visible spectral range: 4.5 - 0.24 pm 

For more information about the RT code go to: http://spacescience.arc.nasa.gov/mars-climate-modeling-group/ 

INITIAL CONDITIONS: 

- Isothermal; 170 K 

- Not connected to 3D model for initial testing 

- Ptrop = 2E-6 mbar; 40 layers 

- Surface pressure = 7 mbar; Albedo = 0.24, Dust free 

- SZA= 80 degrees (for radiative equilibrium case) 


Strategy for ID RT code 

1) Extend tropopause pressure from 8E-4 mbar to 2E-6 mbar (to enter NLTE region) 

2) Increase sigma layers from 24 to 40 (provide more resolution) 

3) Implement NLTE additions 

a) Solar heating rate corrections from Table 1 in Lopez-Valverde et al. (1998) 

b) C0 2 15 pm cooling parameterization adapted from Bougher et al. (2006) 

4) Increase IR spectral bands from 5 bands to 8 bands 

5) Check model to see if it goes to equilibrium and if it is stable 




IR Spectral Bands 


Radiative Equilibrium IR Heat Rates 


1 Band 

Wavelength (pm) 
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1000.0-33.3 
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33.3-18.2 
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18.2-15.7 

Wing of 15 micron C0 2 band 
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15.7-14.2 

Core of 15 micron C0 2 band 
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14.2-12.9 

Wing of 15 micron C0 2 band 
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12.9-10.3 

12 micron water ice band 
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10.3-7.8 

9 micron dust band 
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IR Heating Rates (Individual Bands) 
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All LTE band are represented in panel A (bands 3-5 are scaled by 100 ) 
Bands 3-5 are replaced with the corrected NLTE 15 micron heating rates 
(these bands are shown in panel B) 


Conclusion/Future Work 

Adding NLTE provides smaller heating rates; which is the expected result 

ID RT code does come to equilibrium and it is stable when evolved in time 

More analysis is needed for the upper boundary 

These ID RT code updates will be implemented into the 3D code 

Analysis of MGCM results with the vertical extension and NLTE are still in progress 
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- Output: zonally averaged over 10 sols 
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